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ABSTRACT 

We examine the early angular momentum history of stars in young clusters via 197 photomet- 
ric periods in fields flanking the Orion Nebula Cluster (ONC), 81 photometric periods in NGC 
2264, and 202 measurements of v sin i in the ONC itself. We show that PMS stars spanning 
an age range from 0.1 to 3 Myr do not appear to conserve stellar angular momentum as they 
evolve down their convective tracks, but instead preserve the same range of periods even though 
they have contracted by about a factor of three. This result seems to require a mechanism that 
regulates the angular velocities of young stars. We discuss several candidate mechanisms. The 
most plausible appears to be disk-locking, though most of our stars do not have Iq — K s excesses 
suggestive of disks. However, a decisive test of this hypothesis requires a more sensitive diagnostic 
than the Iq — K s excesses used here. 

Subject headings: stars: pre-main sequence — stars: rotation 



1. Introduction 

Observations over the past decade have established that pre-main sequence (PMS) stars accrete a 
substantial portion of their final mass via material transported through circumstellar disks. However, the 
prediction that accretion of high angular momentum disk material will cause the stars to spin up during the 
accretion phase (e.g. Durisen et al. 1989) is not borne out by observations: most PMS stars have rotational 
velocities (v) of no more than a few tens of km s . By contrast, the typical breakup velocity at ~1 Myr is 
~300 km s" 1 . 

Explanations of the observed slow rotation typically invoke angular momentum transfer either from the 
star to the surrounding accretion disk (e.g. Konigl 1991), or to a stellar wind originating at the boundary 
between the disk and the stellar magnetosphere (e.g. Shu et al. 2000). Either mechanism requires that during 
the disk accretion phase, the angular velocity of the star be "locked" to a period set by the Keplerian angular 
velocity at or near the boundary between the stellar magnetosphere and the accretion disk. At the end of 
the accretion phase, PMS stars should be unlocked from their disks, and free to spin up as they contract 
toward the main sequence. 

Considerable observational effort has been devoted to testing the prediction that rotation periods are 
locked to a narrow range during the accretion phase via comparison of the distribution of stellar periods or 
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rotation rates for samples of stars surrounded by disks with those that lack disks. Early studies suggested that 
stars surrounded by disks tend to have significantly longer rotation periods than their diskless counterparts 
(e.g. Edwards et al. 1993, Choi & Hcrbst 1996); more recent results based on larger samples yield more 
ambiguous results (e.g. Rcbull 2001 [R01]; Stassun et al. 1999 [SMMV]; see also Herbst et al. 2000 [HRHC], 
Carpenter et al. 2001 [CHS], Hcrbst et al. 2001 [HBJM]). HRHC and HBJM both find a bimodal distribution 
of periods for the higher-mass (>O.25M ) stars, whereas it is not bimodal for the lower mass stars. R01 
looked for this dependence, but did not find it. 

Over the past few years, rotation periods P or projected rotational velocities usini have been measured 
for large samples of PMS stars. These data sets are now large enough to to map changes in stellar angular 
momentum (J) as stars of different masses evolve down their convective tracks, and to understand the role 
that disk locking may play in determining the angular momenta of young stars. 

In this paper, we exploit these large datasets to examine the early angular momentum history of PMS 
stars spanning ages from ^0.1 Myr to ^3 Myr and masses ^0.2 to ^2 Mq, both for stars apparently 
surrounded by accretion disks and those that lack evidence of such disks. 

We first examine two samples of low mass PMS stars for which rotation periods are derived from 
observations of spot-modulated light variations: stars located in fields flanking the Orion Nebula Cluster 
and NGC 2264. Together, stars comprising these samples span masses from ^0.2—2 M Q and and nominal 
ages from ^0.1—3 Myr. At fixed T ff, radii appear to span a range of about a factor of 3. Absent any 
external angular momentum loss mechanism (e.g., disk locking or spindown torques exerted by stellar winds), 
contraction of these fully convective stars over this radius range should result in a dramatic (9-fold) decrease 
in period. Hence, our sample appears well-suited to quantifying changes in period distributions among PMS 
stars arising as stars contract along convective tracks, and to assessing the possible role of mechanisms that 
could potentially regulate or mitigate such evolution-driven changes among low mass stars. 

Our approach will be to examine the relationship between observed rotation period (P) and derived 
stellar radius (R) for a large sample of PMS stars - both those which appear surrounded by circumstcllar 
accretion disks (as judged from ncar-IR excess emission), and those that appear to lack such disks. We 
thus first summarize the observational database for each sample: the derived periods as well as sample 
limitations, along with the photometry and spectroscopy that enable locating each star in the Lb \—T c s 
plane (thus yielding a radius) and determining whether or not a star is surrounded by a circumstcllar 
accretion disk. Particular attention is paid to evaluating the uncertainties in derived radii - a critical factor 
in assessing the changes in stellar periods expected as PMS stars contract, absent external regulation. 

We find that PMS stars spanning an age range from ^0.1 to 3 Myr do not appear to conserve stellar 
angular momentum as they evolve down their convective tracks, but instead preserve the same range of 
periods even though they have contracted by nearly a factor of three. Our conclusions apply both to stars 
which appear to lack the near-infrared excesses diagnostic of circumstellar accretion disks and those that 
show conclusive evidence of disks. 

Independent supporting evidence for this result follows both from analysis of the relationship between 
observed periods and derived stellar radii using other recently published surveys, and from recently published 
projected rotational velocities observed for a sample of stars drawn from the Orion Nebula Cluster itself 
(Rhode et al. 2001; RHM). These latter data reveal a decrease of average vsini with decreasing radius, 
consistent with the results derived from the ONC Flanking Fields and NGC 2264 period studies. 

Together, these results seem to require a mechanism that constrains young stars to a constant range 
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of periods from the time they first appear on the stellar birthline to ages at least as great as ^3 Myr. We 
discuss several candidate mechanisms, the most plausible of which appears to be disk-locking. 

2. Observations 

This paper makes use of data sets for three young clusters: the Orion Nebula Cluster (ONC), the 
Flanking Fields (FF) that surround the ONC (cf. Rebull 2001), and NGC 2264. Together, these clusters 
allow us to examine changes in stellar angular momenta for stars with and without IR signatures of disks 
having ages, masses, and radii spanning the ranges <~ 0.1 — 3 Myr, ~ 0.2 — 2 M Q , and ~ 1 — 5 R & respectively. 
In order to place stars in observed and reddening-corrected color-magnitude (CMD) or Hertzsprung- Russell 
(HRD) diagrams, determine stellar radii, and assess whether or not a star has a disk, we require that spectral 
types along with VIqK s photometry be available for all the stars included in the present study. 

We use L bo i and T off to calculate stellar R: 

= (1 > 

To obtain Lbol and T e g from the observed Iq and V — Iq , we first deredden the observed V — Iq using the 
observed spectral type to estimate intrinsic colors. The dereddened Iqo is calculated from the extinction at 
Ic, which is given by Ar=1.61E(V — Ic)- T c g follows from the spectral type. We convert Ico to Lbol using 
the approach described by Hillenbrand (1997; H97); distance moduli were taken from the literature as noted 
in the sections that follow. 

To ensure internal consistency, we recalculated previously published extinction and color excesses for 
each star included in our discussion by combining photometry and spectra reported in the literature with a 
common set of assumed photospheric colors and reddening laws. 

2.1. The ONC 

The ONC is nearby (470 pc; Genzel et al. 1981), compact (size ~1 pc), and young (typical stellar ages 
of ^1—3 Myr; see, e.g. H97). It has been studied recently by Rhode et al. (2001; RHM), who measured 
values of the projected rotational velocity (v sini) for 155 stars (with upper limits for another 83). Observed 
V, Iq, and K s colors and spectral types for these stars are available from H97. There are 202 stars for 
which both the photometry and measurements of i>sini are available; they appear in Figures 1—3. These 
figures show the observed and dereddened CMDs for the stars along with their positions in the L ho \—T cS 
plane (Hertzsprung- Russell Diagram, HRD). About 60% of the stars in the vsmi sample have near-infrared 
(Ic — K s ) excesses indicative of a circumstellar disk, e.g. a de-reddened Ic — K s color exceeding expected 
photospheric values by 0.3 mag, the conservative criterion adopted by H97 (see also Hillenbrand et al. 1998) 
to select the most probable disk candidates. 

A significant advantage of using v sin i data for studies of PMS rotational properties is that line widths 
can be measured for every PMS star to a limit set by spectral resolution; for the RHM data this limit 
is 11 km s" 1 , which corresponds to P > 6 days for a typical PMS star. By contrast, searches for spot- 
modulated periods have biases arising from: (1) the sampling cadence of photometric observations (which 
can be determined if the observation times are known); (2) the greater difficulty in determining periods for 
accreting PMS stars in which stochastic, accretion-driven photometric variability can overwhelm periodic 
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spot-modulated signals; and (3) the possibility that some ranges in stellar rotation rate may not produce 
the large, long-lived spots required to establish periodicity. Together, these biases could well mitigate 
the significant advantage in principle of periods in establishing stellar rotation absent the complication of 
projection effects (sini). 

RHM establish a result essential to studies of stellar angular momenta based on observations of spot- 
modulated periods: there is no statistically significant difference between the vs'mi distributions of PMS 
stars with known periods (P) and those without. Hence there is no statistical bias introduced by determining 
the distribution of rotational properties for that subset of PMS for which period determinations are possible. 

RHM calculate i?sini values by combining observed P and vsinz. They also establish another result 
critical to our analysis: agreement in the mean between these Rsini values and the R values calculated 
from Lboi and T e g. This result supports our use of estimates of R derived from the apparent location 
of an individual star in the HRD as an independent variable that can be used to estimate the effects of 
evolution-driven changes in angular momentum. 

2.2. The Flanking Fields 

The Orion Flanking Fields (FF) were defined and studied by Rebull et al. (2000) and Rebull (2001; 
R01); we use VIcK s photometry, spectral types, and P from these papers, producing 197 stars which appear 
in Figures 1—3. 

Measures of stellar rotation for this sample are provided by spot-modulated photometric periods (P). 
Spot modulation is only found in late-type stars, so unlike vsini measurements, which can be made of 
any type star, P observations are only obtainable for types mid-G and later. Period determinations have 
uncertainties SP <1%. The P measurements used here are complete over a broad range of periods, viz 
P ~ 25 — 0.3d, corresponding to rotational velocities v <~ 3 — 300 km s _1 , ranging from ^1 — 100% of the 
breakup velocity. 

We have limited our primary set of P data in Orion to the studies of the FF cited above in order to 
ensure that we were working with a homogenous data set that was complete over the largest possible range of 
periods. There are several other studies of the Orion region in the literature; they are compared in Table 1. 
Note that the ranges of periods over which these studies are complete (not simply sensitive) can vary widely. 
Results from the other databases spanning more restricted or different ranges of period nevertheless support 
the conclusions reached from analysis of our primary database (R01), as will be shown below. 

Stars with observed periods in the FF are young and thought to be associated with the ONC cluster for 
several reasons. R01 used existing membership studies to demonstrate that the possible field star contami- 
nation in the sample of stars showing periodic variability is low, hence supporting the assumption that stars 
in the FF with derived periods are likely to be associated with the ONC. Further arguments that reinforce 
this contention include that: (1) the CMD of stars in these fields strongly resembles those from the ONC (see 
Rebull et al. 2000 and Figures 1—3), and (2) the spatial distribution of the periodic candidates is elongated 
N-S, along the direction defined by the distribution of dense molecular material which lies behind both the 
ONC and the FF. We thus adopt the same distance, 470 pc, for the FF stars as we did for the ONC sample. 

R01 examines the Iq — K s excesses for the FF sample as part of her study of disk frequency and 
properties. She finds that 30% of the stars in the FF sample have Iq — K s excesses of 0.3 mag (or greater) 
above expected photospheric values - the conservative criterion adopted by H97 and Hillenbrand et al. (1998) 
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as indicative of stars with high likelihood of a circumstellar accretion disk. 



2.3. 



NGC 2264 



Data for NGC 2264 are similar to those used for the FF sample. Colors and spectral types come from 
Rcbull ct al. (2002), and periods from Makidon et al. (2002). A total of 81 stars have spectral types, colors, 
and spot-modulated periods. These stars, appearing in Figures 1—3, come from our study of NGC 2264 
alone; the stars that we have in common with Kearns & Herbst (1997, 1998) are retrieved with identical 
periods. About 20% of the stars in this sample have Iq — K s excesses indicative of a circumstellar disk using 
H97's conservative criterion, Iq — K s excess > 0.3 mag. 

NGC 2264 (part of the Mon OB 1 association) is about twice as far away as Orion; we used a distance 
modulus of 9.40, or 760 pc (Sung et al. 1997). As in Orion, a molecular cloud located behind the cluster aids 
in blocking background field stars and limiting the cluster depth to a small fraction of the distance to the 
cluster. Past studies have argued that NGC 2264 is slightly older than the ONC at 3—5 Myr (e.g. Sung et al. 
1997). Our data suggest that NGC 2264 indeed appears to contain fewer very young stars when compared 
with the ONC, but Figure 2 also suggests that the range of ages found in NGC 2264 is not dramatically 
different from the ONC (Rebull et al. 2002). The value of the NGC 2264 sample for this study thus lies in 
its containing a larger fraction of stars with ages comparable to 3 Myr than either the ONC or FF. 



Key to our study is the ability to estimate stellar R accurately enough to track the changes in P expected 
as PMS stars contract and evolve down their convective tracks. The procedure for calculating R is described 
in the previous section. Sources of error in this calculation include photometric uncertainty, the amplitude 
of stellar variability, uncertainties in the reddening correction, errors in classification, uncertainties in the 
intrinsic photosphcric colors, the effects of accretion, errors in the distance, and the presence of companions. 
All of these effects contribute to errors in Lboi and T e g and hence to uncertainty in R. In the remainder of 
this section, we provide a thorough discussion of uncertainties involved in deriving L ho \, T cff , and R from 
the observed colors and spectral types. Photometric and spectral type errors are summarized in Table 2. 
Unless otherwise specified, Lboi is in ergs cm~ 2 s _1 , T e s in K, and R in R Q . 



3. 



Sources of error in R 



3.1. 



Preliminaries 



Since 



2 log R = log L hol - log(47T(T) - 4 log T eff 



(2) 



and for 



x = au±bv 



(3) 




(4) 



it follows that: 



(<51ogi?) 2 = - A (6\ogL hol f +4(<5 logT cff ) 2 



(5) 
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where we assume that the errors are independent. This assumption gives an upper bound to the total error, 
since errors in spectral type produce partially offsetting errors in T c g and Lboi- 

3.2. Errors in Brightness and Color 

Photometric errors result from photon counting statistics, uncertainties in atmospheric extinction cor- 
rections, and uncertainties in transforming to a standard system. These errors are smallest for the brightest 
stars and largest for the faintest stars, but are typically ^0.03 mag in both Ic and V — Ic for most stars in 
our sample. 

Photometric variability on scales ranging from millimagnitudcs to magnitudes is characteristic of young 
stars (e.g. Herbig 1954). The light curves of the samples with periods (the FF and NGC 2264 stars) have 
typical amplitudes (mean-to-peak) of 0.03 mag in Iq for all but the stars with spectral types M3 and later, 
which have somewhat larger amplitudes. The vsirti sample was specifically selected to include aperiodic as 
well as periodic variables. We do not have monitoring information on these specific stars, but the distribution 
of variations of similar aperiodic stars in the FF and NGC 2264 is strongly peaked near ~0.03 mag as well. 
We therefore take ~0.03 mag as the typical amplitude for all stars in our samples. 

In our photometric studies of the FF and NGC 2264, we monitored in one filter, Ic, only, and thus we 
do not have direct information on how V — Ic or V changes with these periodic modulations. Carpenter 
et al. (2001) monitored stars including those in 2 of the 4 FF in JHK S . They find that most (between 57 and 
77%) of their variable stars can be accounted for by cool spot models, e.g. low-amplitude, nearly colorless 
modulations. Based on this, we assume that the amplitude of variation we found at Iq also occurs in V, and 
that it does not affect V — Ic at significant levels. 

If we assume that V — Ic is affected only by uncertainties resulting from photometric uncertainties, 
for stars with types between K5 and M2, we derive uncertainties in luminosity <51ogLboi ^0.04; similarly 
(51ogT e ff ~ 0.002. Using equation 5 above, SlogR <~ 0.02. 

If we assume that Ic is affected by photometric uncertainties and by variations with an amplitude of 
0.03 mag, we similarly derive, for stars with types between K5 and M2, 51ogLboi ~ 0.02, and SlogT e g <~ 
(because changing Ic alone does not change T e g), and thus SlogR ~ 0.01. These uncertainties increase with 
V — Ic primarily because the redder stars are also fainter. 

3.3. Reddening corrections, photospheric colors, and spectral type uncertainties 

We use spectral types to determine the intrinsic stellar colors and hence to derive (interstellar) reddening. 
The dominant error in this process is the uncertainty in spectral classification. Errors in typing have been 
estimated (H97, Rebull et al. 2000) to be one subclass for the earlier types (F, G, K) and at half a subclass 
for the later types (M). For stars that are still accreting material from their disks, the effects of veiling on 
the observed photospheric spectrum result in classifications that are systematically too early (e.g. H97). 

Another potential source of error is uncertainty in the intrinsic photospheric colors. It is standard 
practice to adopt ZAMS colors as the best estimate for the true photospheric colors for these young, low- 
gravity stars. Our array of ZAMS colors includes measurements from Bessell (1991), Lcggett (1992), and 
Leggett et al. (1998). This issue was discussed in detail by Rebull et al. (2000), so we only summarize some 
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of the most important issues here. (1) Colors for types >M4 are the most uncertain. (2) Chromosphcric 
activity at levels typical of PMS stars can in principle affect PMS colors, but hot active regions are not likely 
to significantly affect colors for passbands V and redder (cf. Rebull et al. 2000). (3) Large dark spots can 
also in principle affect both colors and spectral types; their effects have not yet been quantified. (4) Finally, 
the use of dwarf colors for these young subgiant stars may be questionable. To place an upper bound on this 
latter effect, we have compared colors for dwarf stars with evolved stars of luminosity class I and II. Based on 
this comparison, we note that PMS stars are likely to be slightly redder than dwarfs to about Ml, and bluer 
for stars later than that. However, the size of this error is small compared to the color change effected by 
assuming one spectral (sub-)class uncertainty in each direction. For example, a worst-case scenario occurs 
at M2, where the difference between dwarf and giant V — Iq is ^0.25 mag. By comparison, the difference 
between V — Iq colors of M2 and M2.5 dwarfs is ^0.5 mag. Since the true intrinsic colors of the PMS stars 
are closer to dwarfs than stars of class I and II (cf. H97), this error will be in reality even smaller. 

We have estimated the total error from all of the sources considered in this subsection by shifting the 
observed spectral type one subclass (or half a subclass for the M stars) blue and red, and then comparing the 
resulting Lbol, T e g, an d R. Typical standard errors for stars in the type range K5— M2 are (HogLboi ~ 0.06, 
<51ogT ff ~ 0.01, and 5\ogR <~ 0.04. The uncertainties in T e fj increase as a function of type to mid-K owing 
to the uncertainty in assigning spectral types and the slope of the relationship between spectral type and 
T ff, but decrease from late K through the M- types, because the precision of our assigned spectral types 
improves for these late type stars. The uncertainties in Lboi increase monotonically as a function of type. 

We note that stars in the ONC may suffer from an additional effect: a range in the ratio of total to 
selective extinction that apparently derives from a spatially variable mix of interstellar grain sizes in the 
vicinity of the Trapezium. For a typical V — Iq color excess, the uncertainty driven by a change of a factor 
of two in the total to selective absorption ratio will vary by ~0.5 mag, which corresponds to an uncertainty 
of log iboi of ^0.2 dcx or log R of ~0.1 dex. 

3.4. Accretion 

Material passing through a circumstellar accretion disk is funneled toward the star along magnetosphcric 
columns, and eventually strikes the photosphere at supersonic speeds, producing hot "accretion spots" 
via shock heating. Evidence of these accretion spots is manifest as ultraviolet emission, well in excess of 
photospheric values for stars having large accretion rates. In recent papers, we have exploited excess UV 
emission as one of several measures used to evaluate whether a young star is surrounded by an accretion disk 
(e.g. Rebull et al. 2000, Rebull et al. 2002). For stars with high accretion rates, the accretion spot emission 
can be so strong as to affect V and 7c colors as well. However, stars with accretion spot emission strong 
enough to affect V and Iq are unlikely to be detected as periodic variables owing to the dominance of large, 
random, accretion-driven luminosity changes, as compared with lower amplitude spot-driven variations. 
Quantitative evidence of the effects of accretion on period detectability follows from analysis of the sample 
of pre-main sequence stars in the Orion FF. Of the stars in the present sample with excesses derived directly 
from photometry and spectroscopy, only 30% have detectable periods; by contrast, of the stars likely to 
be cluster members which lack UV excesses, 50% have periods. Consequently, we do not believe that the 
periodic subsamplc considered in this paper is dominated by stars with accretion rates sufficiently high to 
affect V and 7c colors. 

We cannot estimate the role of accretion for the ONC sample. While the majority (60%) of these stars 
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have disks as judged by 7c — excesses >0.3 mag, we have no estimate of the distribution of UV excesses 
among the sample, primarily because the high and variable background in this region of the nebula precludes 
accurate ultraviolet photometry from the ground. 

In this latter context, we note the recent estimates of the effects of accretion-driven changes on PMS 
star colors and derived luminosities carried out by Hartmann (2001). From analysis of a large sample of PMS 
stars in the Taurus- Auriga region, he finds that uncertainties in derived luminosities owing to the effects of 
accretion are typically 0.06 dex in log L bo i, or 0.03 dex in log R. 

3.5. Distance to and width of the clusters 

Uncertainty in distance to the cluster can result in a systematic uncertainty in derived Lboi and R values. 
From the literature, the distance to Orion is 470±70 pc; that uncertainty of 70 pc translates directly to a 
systematic uncertainty of 0.07 dex in log R. Similarly, NGC 2264 is at 760±30 pc, which translates to a 
systematic uncertainty <5 logi?=0.03. Note, however, that an error in the average distance to a cluster simply 
introduces a systematic shift in Lboi or R for all the stars in the cluster. An uncertainty in distance does 
not therefore affect the relative values of R within the cluster, although it will affect comparisons of clusters 
with each other. 

The depths along the line of sight of the three groups of stars studied here appear to be small compared 
to their distances. The apparent angular size of the ONC is (generously) 0.5 deg, corresponding to 1% of 
the distance to the cluster. The scatter in luminosity introduced by variations in the distance to individual 
stars within the ONC is thus much smaller than other sources of error; a 1% error in Lboi corresponds to 
less than 1% error in R. As argued above, in the case of the FF, the presence of the molecular cloud and the 
evident association with and similarity to the young stars in the central ONC all suggest that the distance 
and depth of the cluster in these fields cannot be significantly different from the distance and depth of the 
central ONC. The cluster NGC 2264 also has a molecular cloud behind it and its depth along the line of 
sight is again a small fraction of the distance to the cluster. 

3.6. Binaries 

The presence of binaries could significantly affect the measured Iq and therefore the derived Lboi- The 
worst-case scenario is two equal-mass binaries, where the measured Lboi would be actually twice what it 
should be for a single star; in this case, it would create a 0.3 dex in log Lboi and thus a 0.15 dex change in log 
R. However, the stars studied here are unlikely to have companions with primary to secondary luminosity 
ratios less than about 5:1. For the ONC sample, RHM obtained spectra for all of their stars and would 
have seen double lines either directly or in the cross-correlation peaks had the primary:secondary luminosity 
ratios been smaller than 5:1; they found only 7 binaries in their full sample. For the Orion FF and NGC 
2264, the period sample is also unlikely to be contaminated by unresolved binaries with similar companion 
luminosities. Analysis of such putative binaries would reveal two peaks in the power spectrum analysis of 
their light curves were the primary secondary luminosity ratio smaller than ~5:1. Therefore we estimate 
that the increase in log Lboi as a result of unrecognized binaries is unlikely to exceed 0.10 for any of our 
samples, corresponding to SlogR ~ 0.05. 
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3.7. Final error estimates 

The errors in log R estimated above for each readily quantifiable source are summarized in Table 2 for 
three different mass ranges. In addition, we estimate contributions from accretion (0.03 dex), anomalous 
reddening appropriate to the Trapezium region (0.1 dex), and binarity (0.05 dex). For a typical star in our 
sample (spectral type between K5 and M2), the combination of standard errors from photometry, colors, 
spectra and reddening amount to 0.05 dex (c.f. Table 2). Hence, a reasonable upper limit to the uncertainty 
in log R for stars in this spectral type range would be V0.05 2 + 0.03 2 + 0.1 2 + 0.05 2 = 0.13. 

3.8. Comparison to other error estimates 

These error estimates were derived independently of those discussed in detail by Hartigan, Strom, & 
Strom (1994); our error estimates of <5 log Lboi are identical at 0.08 dex, and our uncertainties in photometry 
and spectral types of 51ogT e ff are essentially identical (0.015 vs. 0.018 dex). 

Hartmann (2001) has also provided a careful assessment of the uncertainties in Lboi and T e g for stars in 
Taurus. His final error estimate for delta log Lboi ranged from 0.09—0.16 dex. However, as noted previously, 
he adopted an error of ^0.06 dex in log Lboi as typical of the uncertainties introduced by the effects of 
accretion on the colors of classical T Tauri studies included in his study. Furthermore, depth for the Taurus 
region compared to its distance is also much larger than for our clusters, adding an additional uncertainty of 
0.06 dex in log Lboi for Hartmann's sample. Even if we were to adopt Hartmann's uncertainties (clearly an 
upper limit to that appropriate for our sample) we would find that the error in log R is less than 0.10 dex. 

3.9. External constraints on errors 

As noted in §2.1, RHM estimate values of stellar radius from (1) log Lboi and T e g (i?(HRD)), and (2) 
from observed periods and projected rotational velocities, vs'mi (i?sini). From 73 stars with both measured 
P and vsini values, we deduce 

logi?sini = 0.93(±0.2)logi?(HRD) - 0.19(±0.1) (6) 

The scatter in log R for each star is ~0.2 dex. Approximately 0.15 dex of this scatter arises naturally from 
the spread in i?sini values introduced by projection effects. The remaining error source must amount to 
~0.13 dex, a value consistent with our upper limit error estimate for stars in the Trapezium region (see 
section 3.7 above). 

Another external check on our error estimates is provided by the analysis of Hartigan et al. (1994). These 
authors estimate luminosities for individual members of PMS binary systems in the Taurus- Auriga region 
using methodology similar to that described above. They then compare the observed ratio L (primary) 
/ L (secondary) with that expected if both stars were born at the same time, and thus fall on the same 
(computed) isochrone. For a sample of 29 stars, two thirds have luminosity ratios within 0.25 dex of 
the expected value, while 90% of the sample lies within 0.4 dex. Hence, subject to the assumption that 
the computed isochrones are correct, this comparison suggests that typical standard uncertainties in log R 
should be 0.12 dex or less. 
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4. Analysis 

If stars conserve (stellar) angular momentum (J) as they evolve down their convective tracks, then we 
would expect that J = Ilu would remain constant (where / is moment of inertia and to is angular velocity) . 
This is equivalent to the requirement that J = MvR be constant if J is conserved in the outermost observable 
layer of the star. Although we expect fully convective stars to rotate as solid bodies, theoretical models of 
PMS stars (e.g. Swenson et al. 1994) predict that changes in surface rotation rate as a star evolves down 
the convective track are the same whether the star rotates as a solid body or conserves J locally (i.e. 
J conservation in spherical shells). The reason is that along the convective track contraction is nearly 
homologous, and changes in I directly track changes in R 2 . For a given mass, if J is conserved, then vR is 
constant and P = 2-kR/v oc R 2 . 

In the following sections, we examine the relationships between R and P and between R and v sin i and 
compare these empirical relationships with the simple prediction that stellar angular momentum is conserved, 
as expected absent angular momentum loss. 

4.1. The FF and NGC 2264 

In Figures 4 and 5, we plot P vs. R for the stars in the FF and NGC 2264. The data are shown for 
three different spectral type groups, which has the advantage of effectively segregating the samples by mass. 
In this way, we can also show the data for stars with spectral types M3 and later separately, since the errors 
in R are largest in this type range. This dividing line is also convenient, as (using D'Antona & Mazzitclli 
models) it corresponds to ~0.25— 0.3 M Q ; this is the location of the break in the distribution of periods as 
found by HRHC and HBJM. 

Reference to Figures 4 and 5 (see also Fig 7) shows that the change in log R is, depending on the cluser, 
at least 0.6 and could be as large as 0.9. Allowing for the standard error of 0.13, it is therefore reasonable 
to assume that there are real changes in log R of at least a factor of 3, a result that is consistent with the 
evidence for contraction obtained by RHM. Figures 4 and 5 show that for all three mass ranges in these 
two clusters, there is no evidence for the spin-up that would be expected with decreasing R and increasing 
age if stellar angular momentum were conserved. Rather, they suggest evolution at constant P, or constant 
angular velocity. 

While Figures 4 and 5 and the fits show no change of average period as PMS stars contract, it is more 
difficult to ascertain via inspection whether the distribution of periods changes significantly. To search for 
such a change, we plot in Figure 6 two histograms depicting the frequency distribution of P for stars in the 
upper and lower quartiles of derived stellar R for our sample. Values of < logP > and < log ii > for the 
upper and lower quartiles are presented in Table 3; note that the average log R for the upper and lower 
quartiles are 0.47±0.01 and 0.13±0.01 respectively, whereas the average log P is 0.50±0.06 and 0.58±0.05. 
Were all stars free to spin up in response to contraction, we would expect the shape of these histograms to 
be the same, but the peak among the lower quartile stars to shift to log P — 0.50 — 0.68 = —0.18. 

If some of the stars were 'locked' to fixed periods, for example by disks, while others were free to spin 
up, we would expect the latter objects to populate an extended 'tail' towards short periods in the solid lines 
in Figure 6. Both visual inspection and quantitative comparison of these distributions via a K-S test reveal 
no evidence that there is a significant difference. The fraction of stars with periods in the lowest (shortest 
period) quartile of the distributions in Figure 6a is 13/52 (25%) for the lower quartile in radius and 19/53 
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(35%) for the upper quartile. Within Poisson statistics, these fractions are identical. 

We conclude from Figure 6 and the above analysis that there is no significant difference in the frequency 
distribution of periods among samples of stars with average log R ~0.47 and 0.13 dex (the upper and lower 
quartile in radius). Hence, there is no evidence of spinup among the bulk of our sample. 

The errors in P are negligible, and so the only way to invalidate this result is to hypothesize that the 
errors in R are much larger than we have estimated-so large, in fact, that we cannot use R to sort stars by 
age. If there were a wide mixture of ages at each R value, then the spin-up expected for conservation of 
angular momentum in diskless stars would be masked by observational scatter. 

Our conclusion that there is no apparent increase in P as stars contract finds additional support from 
analysis of other datasets sensitive to more limited ranges of P and/or R. 

Periods for stars in the ONC can be found in HRHC, SMMV, R01, CHS, and HBJM. We can perform 
a similar analysis of all of these data together for stars just in the ONC, defined for purposes here as the 
region studied by H97. Figure 7 is similar to Figures 4 and 5 but for all 292 stars with available periods, V, 
Ic , and spectral types. Figure 6 presents a similar analysis of the distribution of log P for the upper and 
lower quartiles of log R for this ONC sample; < logP > and < logi? > are noted in Table 3. The fraction 
of stars with periods in the lowest (shortest period) quartile of the distributions in Figure 6a is 21/73 (29%) 
for the lower quartile in radius and 18/73 (25%) for the upper quartile. Within Poisson statistics, these 
fractions are identical to each other and to the results obtained with the FF data. 

We next explore whether these results differ among samples which appear to be surrounded by disks 
and those which appear to lack disks. To enable as robust a comparison as possible, we consider stars from 
the FF sample surrounded by accretion disks to meet the conservative criterion Ic — K s excess > 0.3 mag 
established by H97; those which lack disks are assumed to have Ic — K s excess < 0.1 mag. We exclude 
stars with 0.1 < Ic ~ K s excess < 0.3 mag from the sample on the grounds that establishing the presence 
or absence of a disk is less certain in these cases. Figure 8 is similar to Figure 6 but for stars with Ic — K s 
excesses > 0.3 and < 0.1 mag. As before, Table 3 summarizes < logP > and < logi? >. We conclude that 
there is no significant difference in the distribution of periods among young stars which appear very likely to 
be surrounded by disks and those that lack disks based on Ic — K s excesses alone. This result holds for both 
young stars with large radii and their presumed descendents which have contracted by a factor of three. 

We can fit lines to the log P vs. log R relationships depicted in Figures 4 and 5 (and 7) and to a 
wide variety of subsets thereof. Within the scatter, fits to all of the subsets return a value consistent with 
zero slope. To obtain the largest range of R possible, we combined all of the data from all of the clusters 
and all of the surveys. (The disadvantage to doing this, in addition to the surveys' different sensivities and 
completeness, is that Orion and NGC 2264 are at different distances, and if there are errors in the distance, 
they may affect the comparison.) We broke the entire data set into two samples, with 7c — K s excesses >0.3 
and <0.1. The slope of a fit to these data is 0.2 ± 0.2 for the subsample with no excesses and —0.1 ± 0.2 for 
the subsample with excesses. These results are consistent with a slope of zero (no change of P with R), and 
clearly inconsistent with the slope of 2 expected were stellar angular momentum conserved as stars contract. 

4.2. The ONC 

If stars indeed evolve at constant P to ages t ~ 3 Myr as Figures 4 and 5 suggest, then v/R should be 
constant. In other words, a solid prediction based on the observations of rotation periods in NGC 2264 and 
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the FF is that older stars with smaller radii should be rotating more slowly than younger stars with larger 
radii. If, on the other hand, angular momentum (J) is conserved during evolution, then vR is constant. The 
typical change in R over the populated portion of the convective tracks is about a factor of 3 (cf. Figures 
4 and 5). If we take a typical value of v to be about 35 km s _1 for the youngest stars (see results of RHM 
below), then if J is conserved, the oldest stars should be rotating at more than 100 km s _1 . If angular 
velocity is conserved, then the oldest stars should be rotating at about 12 km s _1 . 

We can make use of the RHM dataset to search for evolution-driven trends in rotational veclocity. To do 
this, we have followed RHM and sorted the stars by Lboi and T c g. Within each (T e ff,Lb i) 'box', we average 
the values of vsini. This approach is necessary because there is substantial scatter in the relationship 
between vsini and R among the stars in this sample arising from (1) the wide range in rotational velocity at 
fixed R (examination of the range of observed periods at fixed radius suggests an expected range of a factor 
of 10), and (2) the range of inclination (sini) values. 

In Figure 9, we superpose the (T e ff,Lboi) 'boxes' on an HR diagram. Also plotted are the stars observed 
by RHM and a set of evolutionary tracks. This diagram shows that the stars observed in the Rhode et al. 
ONC survey span the mass range ^0.2—2.0 M Q . Note also that stars evolve down convective tracks at nearly 
constant T ef j. The RHM approach of sorting the stars by L bo i and L ef j also approximately sorts them by age 
and mass. Each vertical column in the grid shown in Figure 9 spans a range of masses of about 30% around 
the mean mass. The convective tracks are not exactly vertical in this diagram, and so some stars of a given 
mass may move from one column to the next adjacent one as they evolve down their tracks. However, given 
that <vsini> at a given luminosity shows no systematic trend with log T e g (see Figure 10), then the small 
change in average mass with decreasing Lboi within a T e g column will not bias the results significantly. 

As discussed above, the errors in log Lboi and log T e g are estimated to be 0.08 and 0.015 respectively. 
The width of the bins into which the stars were placed are 0.37 in log Lboi and 0.048 in log T e g. It is therefore 
unlikely that the errors in deriving log Lboi and log T e g are large enough to move a star further than to a 
bin just adjacent to its true location. In most of the mass ranges studied, the stars span nearly 1.5 in log 
Lboi, and so the errors in log Lboi are not large enough to mask systematic trends in <vsini> as a function 
of Lboi- Because (cf. Figure 10) there are no systematic trends in vsini with T e g at constant Lboi, errors in 
T e g are not large enough to introduce or mask any systematic trends. 

Figure 10 shows this same rectangular grid, but now with lines of constant radius superposed. Comparing 
this figure with Figure 9 shows that the data for stars in a given mass range typically span a factor of ~3 in 
radius. The average value of vsini for the stars within each of the grid boxes is also shown. The data are 
summarized in Table 4, which gives <vsini>, its standard deviation, and the number of stars in each grid 
box. 

Examination of Figure 10 shows that stars near the tops of the convective tracks have vsini typically 
about 35 km s . Stars near the bottom of the tracks have typical vsini of about 20 km s _1 -/ar smaller than 
the 105 km ,s _1 that would be expected for conservation of stellar angular momentum. Instead, the observed 
value of 20 km s _1 is about 50% higher than the 12 km s _1 that would be expected for conservation of 
angular velocity. However, these results are biased toward minimizing the decrease in vsini with R. In 
calculating these averages, we have used 11 km s _1 for all stars rotating at or below the resolution limit 
set by the spectra analyzed by RHM. Typically 20% or fewer of the stars at the tops of the convective 
tracks are rotating at or below this limit, and so <vsini> for the youngest stars should not be significantly 
overestimated because of the inclusion of stars with measured upper limits. However, more than half 'of the 
stars at the bottoms of the tracks rotate at or below 11 km s _1 , and so the values of <vsini> quoted for the 
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oldest stars are too high. This increase in the fraction of stars with t>sinz< 11 km s is another indicator 
that stellar angular momentum is lost as stars evolve down their convective tracks. 

Note, too, that the decrease in <v sin i> with decreasing radius and increasing age is found independently 
for each of the columns in Figure 10. 

We can estimate the significance of this result by plotting the values of <v sin i> for each of the (T e g,L ho {) 
boxes against the average value < R > for each box; this relationship is shown in Figure 11. Note (c.f. Fig. 10) 
that <vsini> is nearly constant along lines of constant R, and we can therefore combine the data for stars 
of all masses spanned by the RHM sample. From Fig. 11, we derive a formal relationship 

log < vsini >= 0.49(±0.11) log < R > +1.15(±0.05) (7) 

assuming uncertainties of 25% in <vsini> and in log < R >. Upper limits are set to 11 km s _1 in computing 
<t>sini> for each box. 

The positive slope clearly shows that the rotation of stars in this sample slows down as stellar radii 
decrease. The slope of the relation is clearly inconsistent with the slope of —1 that would be expected were 
stars conserving angular momentum as they evolve. 

Table 5 gives for each grid box the product <v sin i> R, which should remain constant down each column 
if angular momentum is conserved, and <wsini> / R, which should remain constant if angular velocity is 
conserved. Within a factor of 2 or better (especially given the 11 km s _1 upper limit on vsini), angular 
velocity is conserved for this ensemble of stars. If we calculate rotational velocities for the FF and 2264 
samples using the derived values of R and observed P, we find the same trend toward decreasing rotation 
with decreasing R as in the ONC. 

Since 60% of the ONC stars appear to have disks (based on having observed Ic — K s excesses > 0.3 mag), 
one might conclude on the basis of this sample alone that we have found strong evidence for disk-locking 
- the most plausible explanation for the apparent decrease in stellar angular momentum with decreasing 
radius/increasing age. As we have seen, however, when we looked at the period data for the FF and NGC 
2264, we found similar evidence for loss of stellar angular momentum in samples of PMS stars that appear 
to be dominated by objects that lack detectable disks using the same detection criteria. 

4.3. Summary of the Observations 

PMS stars appear to preserve the same distribution of periods over radii spanning 1—5 Rq, corresponding 
to ages 0.1-3 Myr. The fact that several different mass groupings of stars in all three clusters show this same 
result, combined with the fact that the independent data sets for the each cluster are characterized by 
different sources of error, appear to make this a robust conclusion. 

If correct, our results require that PMS stars lose angular momentum over this age and radius range. 
We explore several candidate mechanisms below. 
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5. Discussion 
5.1. Disk Locking 

Disk locking provides the most plausible explanation for the above results. However, as noted above, 
we find that the bulk of our sample lack robust evidence of disks, at least as measured by Hillenbrand's 
conservative criterion, 7c — K s excess > 0.3 mag. Smaller values of 7c — K s excess are expected if (1) the 
optically thick accretion disks that produce Iq — K s excesses are viewed nearly equator on; or (2) these disks 
have large inner holes resulting from interruption of the accretion disk by the stellar magnetosphere at radii 
exceeding 5 - a possible outcome if the stellar magnetic field is much larger than the few kilogauss typical 
of TTS or if the disk accretion rate is much smaller than the value of 10~ 8 typical of young TTS. In the 
former case, the radiating area of the disk projected in the observer's direction is small, thus resulting in 
small ncar-IR excesses, while in the latter case, the inner disk lacks dust heated to temperatures sufficient 
to produce a Iq — K s excess. Moreover, observational errors in estimating the 7c — K s excess could result in 
missing a significant fraction of stars with relatively small excesses. 

To explore the possibility that our sample might contain a significant fraction of disks with small near- 
IR excesses, we recalculated the fraction of disks for N2264, the ONC and the FF adopting as a criterion 
7c — K s excess > 0.1 mag. The computed fractions increased from 18, 35, and 58% for 7c — K s excess > 
0.3 to 59, 71, and 78%. This suggests that it is at least plausible that a much larger fraction of our sample 
could be surrounded by disks, and hence that disk locking might well account for the apparent decrease in 
stellar angular momenta as PMS stars contract. 

Support for this speculation is also found from the recent study by Lada et al. (2000). These authors 
undertook an i-band search for disks, noting (1) the much larger excesses above photospheric levels expected 
at 7; and (2) the relative insensitivity of L to the presence of small inner disk holes. They found that about 
45 percent of the M-type stars in the Trapezium had TT-band excesses as judged from their location in the 
J — H/H — K diagram but that about 80 percent of them had excesses at L. 

A second search for disks with large inner holes was conducted by Stassun et al. (2001), who searched 
for 10 /xm excesses in stars in Taurus- Auriga and Orion that lack 7c — K s excesses. In this case, only 3 of 32 
stars showed evidence of 10 /im excesses, and on this basis the authors argue that most of the stars without 
near-IR excesses do indeed lack disks. 

From these two studies, it is therefore difficult to estimate what fraction of stars surrounded by disks 
are missed using 7c — K s excess as a disk discriminant. 

Thus, while it appears tempting to argue that disk-locking accounts for the apparent immutability of 
the period distribution among contracting PMS stars, more work is needed in order to establish samples of 
PMS stars for which periods are available and for which there is definitive evidence establishing the presence 
or absence of accretion disks. 

SIRTF observations promise the sensitivity for rapid and accurate 3—10 /xm surveys capable of identi- 
fying disk candidates unambiguously. SIRTF should also establish the (unlikely but still plausible) existence 
of gaseous accretion disks in which all but a fraction of an Earth mass of the small dust grains responsible 
for producing excess IR emission have been agglomerated into larger bodies. 
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5.2. Disk-Locking Combined with Birthline Effects 

Could the combination of similar period distribution for stars of apparently different ages/radii and the 
lack of clear distinction in the period distribution among stars with and without disks be explained through 
some other mechanism? 

As one possibility, suppose that (1) stars in Orion and NGC 2264 were born in a single burst of star 
formation ~1 Myr ago; (2) the range in L (equivalently R) for the PMS stars in these clusters docs not 
reflect evolution down convective tracks, but rather differences in mass accretion rate (M) that force stars to 
evolve along different 'birthlincs' and therefore to arrive at different initial luminosities along the convective 
tracks (Palla & Stahler 1992); and (3) PMS stars are locked to a particular P so long as they are surrounded 
by accretion disks. In this picture, the observed distribution of stars along a convective track for a given 
mass then reflects a range of M: protostars with higher M have larger initial radii and lie higher in the 
color-magnitude diagram following the end of the accretion phase (see also Hartmann et al. 1997). The 
similarity of P at different R then follows from the assumption that stars are locked to their disks for much 
(nearly all) of their ~1 Myr accretion phase. 

The main difficulty with this proposal lies with the additional requirement that stellar L cannot have 
decreased much since the stars were deposited on their convective tracks. If stars that formed through high M 
and were initially deposited high on their convective tracks had subsequently evolved downwards, we would 
expect to see a 'tail' of stars that have spun up to shorter P mixed in with the ensemble of objects that 
started their evolution at smaller R. We do not (see Figure 6). The only way to resolve this contradiction 
and still maintain conservation of stellar J would be to identify a mechanism for halting or slowing evolution 
of stars deposited high on their convective tracks for a time comparable to at least ~3 Myr. 

5.3. Angular Momentum Loss via Stellar Winds 

Stellar winds loaded onto open magnetic field lines can exert a spindown torque on stars. However, this 
mechanism appears to be ineffective for PMS stars because, during this phase of evolution, the timescale for 
spindown exceeds the evolutionary timescale by a few orders of magnitude (e.g. MacGregor and Charbonncau 
1994) - fully convective stars are assumed to rotate as solid bodies, and the wind must slow down the entire 
star. 

Calculations (e.g. Kawaler 1988) show that the rate of change of J depends on the product of the mass 
loss rate and the square of the Alfvcn radius, but the Alfven radius varies inversely as some power of the 
mass loss rate, with the specific power depending on the configuration of the magnetic field. For a field 
geometry "intermediate" between a dipolar and a radial field, dJ/dt does not depend on mass loss rate at all 
(Bouvier et al. 1997). The only circumstance under which magnetic winds could play an important role in 
slowing the rotation of PMS stars would be if there were some way in which to decouple the outer layers of 
the star from the interior, i.e. if J were conserved locally, and the wind had to slow down only a thin outer 
layer of the star. Such decoupling is not expected for fully convective stars. 

A further challenge to any wind-driven J loss mechanism is the additional requirement that significant 
J loss would have to cease on timescales of no more than a few Myr in order to account for the significant 
population of stars that ultimately arrive on the main sequence as the rapid rotators observed in young 
clusters such as aPersei (e.g. Stauffer et al. 1989); such stars require spinup from the PMS to the ZAMS. 



- 16 - 



5.4. Angular Momentum Loss via Tidally-Locked Planetary Companions 

Recent theoretical models indicate that tidal locking between a close-in Jupiter-mass planet and the 
parent star can transfer spin J from the star to the orbital J of the planet, thus slowing the rotation of the 
star while driving the planet into a larger orbit. 

This effect is not large enough to account for the observations reported here. The difficulty lies in the 
apparent incompatibility of two requirements: (1) the putative planet must be located close enough to the 
star to produce significant tidal distortion; and (2) the planet must be located far enough from the star to 
dominate the J of the system (and thus be able to create significant changes in stellar J with only modest 
orbital evolution). A simple calculation (cf. Trilling et al. 1998) in which we assume a Jupiter-mass planet 
orbiting a PMS star at a distance of 0.1 AU suggests that J regulation by tidal locking to such a planet fails 
by 2 orders of magnitude. 

6. Summary 

The observations of PMS stars discussed here span an age range of about 3 Myr and a change in stellar 
radius of about a factor of 3. Surprisingly, observations of both rotation periods and projected rotational 
velocities for a sample of several hundred stars show that most of the stars, whether surrounded by observable 
disks or not, preserve constant angular velocity ui to within a factor of 2 as they evolve downward along their 
convective tracks - and thus must lose stellar angular momentum J in direct proportion to R 2 . 

A recent theoretical paper by Tinker, Pinsonneault, & Terndrup (2002) also concluded that stars must 
lose significant amounts of angular momentum during the first few million years after they arrive at the 
birthline. These authors assumed that the ONC stars, when they reach an age of 120 Myr, should produce 
a rotational velocity distribution that looks like the distribution in the Pleiades, and then calculated how 
the angular momentum of the ONC stars would have to change with time to achieve this outcome. Tinker 
ct al. concluded that there must be a mechanism in addition to stellar winds that produces loss of stellar 
angular momentum. They further showed that, if this mechanism is disk-locking, then it must be effective 
for all stars up to an age of about 3 Myr, or alternatively disks must have a range of lifetimes over which 
they serve as effective brakes, and that these lifetimes extend up to 6 Myr in some stars. 

These theoretical ideas combined with recent estimates for the lifetimes of disks may explain why we 
have not yet detected evidence for spin up in the samples we have observed. In particular, Haisch, Lada, 
& Lada (2001) show that more than 80% of the stars in NGC 2024, which is about 10 6 years old have IR 
excesses indicative of disks, while only 50% of the stars in NGC 2264, which they estimate to have an mean 
age of 3xl0 6 years, have such excesses. Assume that all stars with ages 10 6 years and less are surrounded 
by disks and that during this time, their rotation periods are locked. Suppose further that half of these stars 
lose their disks sometime between 10 6 and 3xl0 6 years. Reference to Figures 4, 5, and 7 shows that fully 
half of our sample stars have ages ~1 Myr or less and thus are, by hypothesis, locked to their disks. The 
remaining stars in our sample have ages between ~1 and ~3 Myr. Of those stars, at least half should still 
be surrounded by disks. Hence less than a quarter of our sample might be free to spin up, and most of these 
objects are likely to have lost their disks within the past 1 Myr, and not have spun up significantly. If these 
estimates of disk lifetimes are correct, it seems plausible that for the distribution of ages represented among 
our sample stars and the likely distribution of disk lifetimes, that such modest spinup might well be masked 
by the intrinsic spread in P characteristic of the youngest (presumably disk-locked) stars in our sample. 
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The essential point, however, is that both the observations reported here and theory seem to require 
that a braking mechanism-whatever it may be-must be effective for most pre-main sequence stars during at 
least the first 3 Myr or so after they reach the birthline. 

Significant spin-up must ultimately occur in at least some stars in order to account for the rapid rotators 
seen on the main sequence in young clusters such as Alpha Persei (Stauffer et al. 1989). However, recent 
observations of vs'mi for a small sample of stars in the 10 Myr old TW Hya association (Torres et al. 2001; 
Sterzik et al. 1999) suggest mean vs'mi values similar to those found among the oldest stars in our NGC 
2264 sample, despite the fact that the nominal radii among the TW Hya stars are smaller by 50%. These 
observations suggest that "w locking" may extend to 10 Myr and makes urgent the need for a campaign 
focused on mapping the angular momentum evolution of low mass stars in the age range 3—30 Myr so as to 
understand when rotational spinup in response to contraction takes place. Additionally, it would help to have 
an age indicator other than radius or else to show that radius correlates at least roughly with independent 
age estimates for clusters other than the ONC. 
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Table 1 : Surveys of stars in the Orion region 



description 



SMMV 



HRHC a R01 a 



CHS a HBJM a 



# stars monitored 
monitoring band 

# periods reported 
range of sensitivity (d) 
range of completeness (d) 
approx. mass range (Mq) 

# stars used here d 



2279 

Ic 
254 
0.1-10 
0.1-6 
0.15-1.0 
142 (93) 



134 281 
l-35? c 0.2-40 
l-? c 0.3-20 

0.2-2.0 0.2-1.3 

95 (30) 220 (173) 



829 3585 
Ic Ic 



17,808 2294 

JHK I c 

233 369 b 

1- 13 0.2-30? c 

2- 10 0.8-15? c 
0.2-1.0 0.1-1 
70 (16) 211 (145) 



a SMMV=Stassun ct al. (1999), HRHC=Herbst et al. (2000), R01=Rebull (2001), CHS=Carpenter et al. (2001), HBHM = Hcrbst 
et al. (2001) 

6 404 stars were reported in HBJM, but 369 remain in the most current list provided by W. Herbst. 

c The range over which this study is complete is not reported, and completeness cannot be determined without specific times of 
observations, which are also not reported. 

d The most significant limitation on the databases is that we require that spectral types and VIcK a photometry be available 
for all the stars included in the present study. However, because of significant overlap in each of these Orion fields, there are 
of course duplicate stars in these databases. The second number in parentheses in each column denotes the number of stars 
actually used here from the corresponding database. For example, there are 142 stars in the SMMV database for which SMMV 
measured a period, but only 93 for which they alone measured a period (or have the best period available). In the remaining 
49 stars, another more recent study retrieved an identical period, or has made a convincing case that it has the better period 
for those stars. 
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Table 2. Mean magnitude of empirically-derived errors 



category 


uncertainty 


<K5 


K5— M2 


>M2 


photometric 


Sloe T e ff 


0.004 


0.002 


0.002 


in V - I 


(5l0£ Z/K^l 
w ivy O J - J D01 


0.052 


0.044 


0.039 


(§3-2) 


<51og i? 


0.027 


0.022 


0.020 


photometric 


<51og T eff 











in I + 


(Slog L bo i 


0.017 


0.025 


0.022 


amplitude (§3.2) 


SlogR 


0.009 


0.013 


0.011 


spec type 


Slog T eff 


0.013 


0.012 


0.010 


(§3-3) 


(51og L b oi 


0.047 


0.063 


0.119 




Slog R 


0.035 


0.040 


0.063 


net mean 


Slog T e ff 


0.014 


0.012 


0.010 




Slog Lboi 


0.072 


0.081 


0.127 




<51og i? 


0.050 


0.050 


0.070 


typical 


<51og Teff 




0.01 






(nog L b oi 




0.08 






<51og i? 




0.05 





Table 3. Mean log P and log R for upper and lower quartiles in log R (see text) 



sample 


upper < logP > 


upper < log R > 


lower < logP > 


lower < log R > 


KS prob a 


X 2 prob a 


FF 


0.50±0.06 


0.47±0.01 


0.58±0.05 


0.13±0.01 


0.22 


0.74 


ONC 


0.61±0.04 


0.54±0.01 


0.57±0.04 


0.15±0.01 


0.54 


0.63 


FF disks b 


0.47±0.09 


0.53±0.02 


0.62±0.09 


0.16±0.03 


0.07 


0.20 


FF no disks 


0.58±0.10 


0.41±0.02 


0.46±0.09 


0.10±0.02 


0.30 


0.37 



a The probability derived from a two-sided K-S test of the continuous (or x 2 test of the binned) distributions 
of log P for the upper quartile of log R compared with that for the lower quartile of log R. Neither the x 2 nor 
the KS test suggest any statistically significant difference between the upper and lower quartiles. 

b Those stars with Iq — K s excesses > 0.3 mag 

c Those stars with Ic — K s excesses<0.1 mag 
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Table 4. Distribution of vsmi for the ONC stars. 



log T cff a ^ 


A (3.730-3.682) 


B (3.682-3.634) 


C (3.634-3.586) 


D (3.586-3.538) 


E (3.538-3.490) 


log L/LqI 


<vsmi> 


num 


<vsmi> 


num 


<vsmi> 


num 


<v sinz> 


num 


<vsini> num 


1 (1.48 - 1.11) 


33±13 


3 


41±10 


5 












2 (1.11 - 0.74) 


30±4 


6 


26±6 


4 


14 


1 








3 (0.74 - 0.37) 


26±5 


7 


15±4 


4 


34±12 


7 


35±9 


5 




4 (0.37 - 0.00) 


19±7 


2 


17±3 


13 


17±3 


18 


23±4 


17 


34±15 3 


5 (0.00 - -0.37) 


12 


1 






20±8 


6 


19±3 


41 


19±2 22 


6 (-0.37 - -0.74) 










18±7 


3 


16±4 


18 


20±3 44 



a Subdivisions defined exactly as in RHM; subdivisions according to T e g and L correspond to subdivisions by M for 
stars on convective tracks, and with age for stars of a given M. Column A from RHM spans a larger range of masses 
than the other columns, which span log M/M & ^0.1 dex when calculated via D'Antona & Mazzitelli (1994) or Siess 
et al. (2000). 

b Number of stars in the bin. 



Table 5. Trends in angular momentum (cx vR) and angular velocity (oc v/R) for the ONC. 



log T cff a ^ 


A (3.730-3.682) 


B (3.682-3.634) 


C (3.634-3.586) 


log L/LqX 


< vsmi > /R° 


< v sin i > x R° 


< v sin i > /R < v sin i > x R 


< vsmi > /R 


< v sin i > x R 


1 (1.48 - 1.11) 


6 


191 


6 296 






2 (1.11 - 0.74) 


8 


112 


6 123 


2 


82 


3 (0.74 - 0.37) 


11 


64 


5 46 


9 


131 


4 (0.37 - 0.00) 


12 


31 


8 34 


7 


43 


5 (0.00 - -0.37) 


11 


13 




12 


33 


6 (-0.37 - -0.74) 








17 


19 



Subdivisions defined exactly as in RHM; subdivisions according to T e g and L correspond to subdivisions by M for stars < 
A from RHM spans a larger range of masses than the other columns, which span log M/M ~0.1 dex when calculated via 1 

b < vsmi > /R would be constant if stellar J is conserved, and < vsmi > /R would be constant if lo conserved. Note 
within a factor of <~2. 
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Fig. 1. — Observed color-magnitude diagram (CMD) for each of the three samples considered in this paper: 
the Orion Nebula Cluster itself, the Orion Nebula Cluster Flanking Fields, and NGC 2264. Iq magnitudes 
have been converted to absolute magnitudes (see text for adopted distance moduli) for ease of comparison 
between the Orion samples and NGC 2264. ZAMS relation is heavy line, and V — Iq colors corresponding to 
various spectral types are indicated in top panel. In bottom panel, typical standard error bars are indicated 
for points with V — Iq between 0—1, 1—2, 2—3, & 3—4; see text. 
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Fig. 2. — Dereddened CMD for each of the three samples considered in this paper: the Orion Nebula Cluster 
Flanking Fields, the Orion Nebula Cluster itself, and NGC 2264. Iq magnitudes have been corrected to 
absolute magnitudes for ease of comparison between the Orion samples and NGC 2264. V — Iq values 
appear "quantized" because stars have been reddening-corrected to the intrinsic V — Iq colors expected for 
their spectral types. ZAMS relation is heavy line, and V — Iq colors of spectral types are indicated in top 
panel. In bottom panel, typical standard error bars are indicated for points with V — Iq between 0—1, 1—2, 
2—3, & 3—4; see text. 



-24- 




Fig. 3. — Theoretical Hcrtzprung- Russell Diagrams (HRD) for each of the three samples considered in this 
paper: the Orion Nebula Cluster Flanking Fields, the Orion Nebula Cluster itself, and NGC 2264. Lboi and 
T e ff have been derived from dereddened colors following the formulas given in Hillenbrand (1997). ZAMS 
relation is heavy line, and V — Iq colors corresponding to various spectral types are indicated on top panel. 
In bottom panel, typical standard error bars are indicated for points with log T e g between 3.5—3.6, 3.6—3.7, 
& 3.7-3.8; see text. 
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log (R/Fy 

Fig. 4. — Period vs. radius for stars in the FF. (a) stars with types K5 and earlier; (b) stars with types 
K5-M2; (c) stars with types M3 and later. Open circles are stars without Ic — K s excesses (i.e. without disks) 
and solid circles are stars with excesses (i.e. with disks). Typical standard error bars for each type range are 
indicated. The line in each panel indicates the slope expected if stellar angular momentum is conserved. 
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Fig. 5. — Period vs. radius for stars in NGC 2264. (a) stars with types K5 and earlier; (b) stars with types 
K5-M2; (c) stars with types M3 and later. Open circles are stars without Ic — K s excesses (i.e. without disks) 
and solid circles are stars with excesses (i.e. with disks). Typical standard error bars for each type range are 
indicated. The line in each panel indicates the slope expected if stellar angular momentum is conserved. 
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Fig. 6. — Histograms of distributions of periods for stars in (a) the FF and (b) the ONC for the upper and 
lower quartiles in radius. There is no difference in the period distributions. See Table 3 for < logP >, and 
< logi? >. 
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Fig. 7. — Period vs. radius for stars from the ONC, defined as the region covered by H97, periods taken from 
HRHC, SMMV, R01. CHS, and HBJM. (a) stars with types K5 and earlier; (b) stars with types K5-M2; 
(c) stars with types M3 and later. Open circles are stars without Iq — K s excesses (i.e. without disks) and 
solid circles are stars with excesses (i.e. with disks). Typical standard error bars for each type range are 
indicated. The line in each panel indicates the slope expected if stellar angular momentum is conserved. 
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Fig. 8. — Histograms of distributions of periods for stars in the FF with (a) Ic — K s excess >0.3 mag and 
(b) <0.1 mag. There is no difference in the period distributions. See Table 3 for < logP >, and < logi? >. 
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Fig. 9. — The HRD for the ONC stars. The evolutionary tracks are from D'Antona and Mazzitelli (1994). 
The grid shows the bins in the log Lboi and log T e g plane into which we sorted the stars in order to look 
for trends in <t>sini> with decreasing luminosity and increasing age. The estimated error bars for log Lbol 
and log T e ff are given in the lower left corner of the plot. 
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Fig. 10. — The values of <vsini> in km s -1 as a function of position in the log Lboi and log T e g plane. 
Lines of constant radius are shown as dashed lines. Note that the trend in all columns is toward decreasing 
<wsini> with decreasing luminosity and decreasing radius, which indicates that angular momentum must 
be lost as stars evolve down their convective tracks. Table 4 indicates how many stars fall within each box 
in the grid and gives the statistical errors in the values of <wsini>. Note that there is only one star in the 
top box in the middle column. 
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Fig. 11. — Mean values of vsini and log R for each of the (T e ff,iboi) 'boxes' from the previous plots. Note 
from previous figures that <vsini> is nearly constant along lines of constant R, thus encouraging us to 
combine the data for stars of all masses spanned by the RHM sample. The linear fit indicated here is 
log < vsini >= 0.49(±0.11) log < R > +1.15(±0.05). 



